Entanglement of Formation for Gaussian States 
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The entanglement of formation (EOF) is computed for arbitrary two-mode Gaussian states. Apart 
from a conjecture, our analysis rests on two main ingredients. The first is a four-parameter canonical 
form we develop for the covariance matrix, one of these parameters acting as a measure of EOF, 
and the second is a generalisation of the EPR correlation, used in the work of Giedke et al [Phys. 
Rev. Lett. 91, 107901 (2003)], to noncommuting variables. The conjecture itself is in respect of an 
extremal property of this generalized EPR correlation. 
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Entanglement is an essential resource for many quan- 
tum information processing tasks, and hence it is impor- 
tant to be able to quantify this resource. A reasonable 
set of demands lead, in the case of bipartite pure states, 
to a simple and unique measure for this resource: it is 
the von Neumann entropy of either subsystem ^li [3j ■ 
For mixed states, however, many different entanglement 
measures continue to be under consideration [3| . One of 
these measures with an attractive physical motivation 
is the entanglement of formation (EOF) Q. The asymp- 
totic version of EOF is the entanglement cost \3, ^ . EOF 
is defined as an infimum: 

EOF(p)^inf{^p,i?(V',) I P = Ep^IV'.)(^^I}- 

j j 

The infimum is to be taken over all realizations of the 
given mixed state p as convex sum of pure states, and 
E{ipj) = S{tTB[\'^Pj){'4'j\]) where S{-) is the von Neumann 
entropy. EOF has been computed in closed form for arbi- 
trary two-qubit states [Tl), and for highly symmetric states 
like the isotropic states |q| and the Werner states Q. 

Gaussian states, whose entanglement originates in non- 
classicality of the squeezing type , have played a dis- 
tinguished role in quantum information in respect of 
continuous variable systems 
tation 



Their use in telepor- 



12, Ui 



and quantum cryptography [lj| has been 
demonstrated. Q uestions related to their separabil- 
ity 15, [3, 17, [13I and distillabihty [IB] have been re- 
solved. More recently, analytic expression for their EOF 
has been obtained in the symmetric case (20I] . This no- 
table achievement seems to be the first computation of 
EOF for states of infinite rank. These authors exploit a 
certain extremality that the two-mode-squeezed vacuum 
enjoys in respect of the Einstein-Podolsky-Rosen (EPR) 
correlation [2l| on the one hand and entanglement on the 
other. Further analysis of EOF in this case has been 
made j2^] from the viewpoint of Bures distance 

An interesting Gaussian-state-specific generalisation of 
EOF, the Gaussian_entanglement of formation, has also 
been explored i23,i24| . But the EOF of asymmetric Gaus- 
sian states has remained an open problem jl^ . 

In this Letter we compute, under a conjecture, the 



EOF for arbitrary two-mode Gaussian states. Our anal- 
ysis rests on two principal ingredients. The first one is 
a four-parameter canonical form we develop for the co- 
variance matrix; one of these parameters proves to be a 
measure of EOF. The second one is a family of gener- 
alised EPR correlations for noncommuting pairs of non- 
local variables; this family is indexed by a continuous 
parameter 0. And the conjecture is in respect of an ex- 
tremal property of this generalised EPR correlation. 
Canonical Form for Covariance Matrix: Given a two- 
mode Gaussian state, with the mode on Alice's side de- 
scribed by canonical quadrature variables xa, pa and 
that on Bob's side by xb, Pb, we can assume without 
loss of generality that the first moments of all four vari- 
ables vanish 1^, 2^] . Such a zero- mean Gaussian state is 
fully described by the covariance matrix 3) 23] 
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where the phase space variables are assumed to be ar- 
ranged in the order {xatPAtXb,Pb) = S,, and we have 
retained through the parameters a, [3 > Q the freedom 
of independent local unitary (i.e., symplectic) scalings on 
the A and B sides. This freedom will be used shortly. 

Note that Vg is left with no correlation between the 
'spatial' variables xa, xb and the 'momentum' variables 
PA, Pb- Thus it is sometimes convenient to view Vg as 
the direct sum of 2 x 2 matrices: 



Vg=Xg®Pg, 



Xg = ^ 



an 



a 
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Let I'^r) denote the standard two-mode-squeezed vac- 
uum state with squeeze parameter r. R takes the 
Schmidt form in the standard Fock basis: 



c„]n)^ (g) \n 

n=0 



/ ^ ^n\'-/ A I"'/ b — ^ ^ c„]n, n) , 

n=0 n=0 

tanh"r/ coshr . (2) 
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Denoting by Er the entanglement of \'^r), we have 

Er — cosher log2(cosh^r) — sinh^r log2(sinh^r) . (3) 
The covariance matrix of |^', ) has the form 
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■ C 


s' 




' c -s' 
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' p-^ = \ 


-s c 



C = cosh 2r, S = sinh 2r . 



(4) 



Proposition 1 : Given a two-mode covariance matrix Vg, 
the local scale parameters a, (3 can be so chosen that Vq 
gets recast in the form 



C + uc^ S + ucs 

C + vc^ -S'-i;cs 

5' + ucs C + us^ 

-S'-vcs C + vs^ 

C = cosh2ro, S' = sinh2ro; c = cos0Oj s = sin0o- 



Note: We will call Vq the canonical form of a two-mode 
covariance matrix; our results below will justify this el- 
evated status. We assume without loss of generality 
n > m or, equivalently, < 9o < 7r/4. For a given Vg 
there will be two solutions for the above form. Canonical 
form will always refer to the one with the smaller squeeze 
parameter ro, which is ensured by the restriction 



tan 6q > tanh rg . 



(5) 



This condition proves central to our analysis. Its origin 
may be appreciated by inverse two-mode-squeezing the 
Gaussian state Vq until it becomes just separable, and 
noting that there exists a range of further squeezing in 
which the mixed Gaussian state remains separable before 
becoming inseparable again. The parameters m, u > 0. 
The essence of the canonical form is that Vq differs from 
the covariance matrix of a two-mode- squeezed vacuum 
I'^ro) by a positive matrix which is a direct sum of two 
singular 2x2 matrices which are, modulo signature of 
the off- diagonal elements, multiples of one another. 
Proof: The canonical form demands, as a necessary con- 
dition, that a, (3, and r be chosen to meet 

Aet{XG - X^^) = , dct(PG - f*,.) = . (6) 

These being two constraints on three parameters, one 
will expect to get a one-parameter family of solutions to 
these constraints. For each such solution we may denote 
the vector annihilated by the singular matrix Xq — X^,^ 
by (sin 6*, — cos0), and that annihilated by Pq — P*,- by 
(sin 6* , cos 6* ). The canonical form corresponds to that 
solution for which = 0; it is this degenerate value that 
equals Oq of the canonical form. 

That there exists such a degenerate value can be seen 
as follows. We may fix the scale parameter a through 



a = ^Jmjn, and then solve Eqs. (6) for (3 and r, the 
smaller r being the relevant one. We will find 9 = -k/A 
and 9 < 7r/4 in this case. On the other hand if we take 
a = ^Jn/m and then solve Eqs. (6), we will find 9 — 
7r/4 and 9 < it/ A. It follows from continuity that there 
exists an intermediate value ao for the parameter a, in 
the range ^Jm/n < a < ^Jn/m, for which 9 = 9 {< n/A 
since n > to). And this yields the canonical form. 

Viewed alternatively, the canonical form Vq places the 
following two requirements on the scale factors a, (3: 



detXG-1/4 tT:{a-iXG) 



detPc - 1/4 tr(a3PG) 

det(XG-cr3PGf73) = 0, 



(7) 



where 173 is the diagonal Pauli matrix. These are simul- 
taneous equations in a, (3, and solving these equations 
yields, in terms of n, m, kx, kp, the values of a, (3 corre- 
sponding to the canonical form. 

Two special cases may be noted. If to = n we 
have a = \ (since ^Jn/m = ^Jm/n), and hence (3 = 

[n — kp)/{n — kx), so that the canonical squeeze pa- 
rameter To is given by e"^''" = {n — kx){n — kp), re- 
producing the results of Ref. [20] . The parameter 6*0 al- 
ways equals 7r/4 in this (symmetric) case. On the other 
hand, li kx — kp — k, the canonical form corresponds to 
a = /3 = 1, and one obtains by simply solving 



det 



n — cosh2ro k — sinh2ro 
k — sinh2ro to — cosh2ro 



0, 



(8) 



which yields this closed- form expression for ro: 

nra — k^ + 1 



cosh(277 — 2ro) 



.±2r, 



y^{n + niY - Ak"^) 

(n + to) ± 2k 
a/(to + n)2 - 4fc2 



Generalised EPR Correlation : To proceed further, we 
need to generalise the familiar EPR correlation [2^ . 
Given any bipartite state jV'), define 



xg = sm9 XA — cos9 XB, Pe - 



sm9pA -\- cos9pB , 

^mpom. (9) 



In defining Ag{ip) we have assumed {ip\xe\^) = = 
{iplpelip}) if this is not the case then xq and pg in Agiijj) 
should be replaced by xg — ('0|xe|-0) and pg — MpoM 
respectively. Clearly, the usual EPR correlation [20| cor- 
responds to 9 ^ t:/A. While Xy^/^, PT^/i commute, the 
generalised EPR (nonlocal) variables xg, pg do not com- 
mute, and hence the name generalised EPR correlation 
for Kg{'^); indeed, we have [xg, pg] — —icos29. For the 
two- mode-squeezed vacuum j^'^) the generalised EPR 
correlation reads 



Ag{'^r) = cosh2r - sin 26^ sinh 2r . 



(10) 
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Let us combine the quadrature variables of the os- 
cillators of Alice and Bob into boson operators a = 
{xA + ipAj/V^ and b = (XB +ipB)/V2. Then, Ae(V') 
has this expression quadratic in the boson variables: 

Ae(v) = imow, 

-2cosesm9{ab + a''b^). (11) 

We may call Ag the generalised EPR operator. 

The entanglement of \^r) monotonically increases with 
increasing value of the squeezing parameter r. In order 
that Ae(^'r) be useful as an entanglement measure of \^r) 
it should, for fixed value of 9, decrease with increasing 
r. The restriction tan0 > tanhr, encountered earlier 
in Eq. (5) from a different perspective, simply ensures 
this. Through the monotonic relationship (3) between 
r and Er, we will view this constraint as a restriction 
on the allowed range of values of 6', for a fixed value of 
entanglement. 

Given a squeezed state |5'r), let us denote by \^'^) the 
state obtained from \^r) by independent local canonical 
transformations fl^ Sa^ Sb € Sp{2,R), acting respec- 
tively on the oscillators of Alice and Bob. 
Proposition 2: We have Ae(^'J.) > Ae(*r), in the 
range 1 > tan 6* > tanhr and for all Sa, Sb & Sp{2, R). 
Proof: Clearly, Ae(1';) = \ {cosk2r[sm^ 9 tv{SAS'^) + 
cos2 9 iT{SBSl) ] - sin 261 sinh 2r tr (aa^^CTsS'l) }. If e^^"^ 
are the singular values of Sa, and e^^^ those of Sb, 
then tr(5'AS'J) = 2cosh27A, tr(S'BS'|J) = 2cosh27B, 
and tT:{(JzSA<^'iS]^) < 2cosh(7A +7b). Thus the differ- 
ence A(7A, 7s) = Ae{%) - Aei'fr) obeys A{-fA, 7b) > 
cosh2r[sin^ 6'(cosh27A - 1) + cos^ 6l(cosh27B - 1)] - 
sin26'sinh2r[cosh(7A -I- 75) — 1]. It is easily seen that 
A(7A, 7b) is extremal at 7^1 = 75 = corresponding to 
the standard squeezed state |^r). To show that this ex- 
tremum is indeed minimum we note that the determinant 
of the Hessian matrix of the right hand side, evaluated at 
7a = = 7b, is proportional to sin 20 cosh 2r — sinh2r, 
and hence is positive if and only if tan^ > tanhr. 

Once again we see a role for the requirement tan 9 > 
tanhr. Let the equivalence ~ Vq denote the fact 
that the corresponding Gaussian states are connected by 
a local canonical transformation. The fact that M = 
Vo - V*^^ > implies Aeg{pVo) > ^Ooi'^ro)- In view 
of Proposition 2 this implies Ag„{pva) > Ag„{pva) > 
^doi'^ro) = cosh2ro — sin 26* sinh2ro for any Gaussian 
state Vg connected to Vo by local canonical transforma- 
tion. This assigns an alternative meaning to the canoni- 
cal parameter rp: 

Proposition 3 : Given a Gaussian state described by Vg ^ 
Vo, the canonical squeeze parameter rg is the smallest r 
for which the matrix inequality Vg — V^i > is true. 

It is well known that the two-mode-squeezed vacuum 
has several extremal properties of interest to entangle- 
ment [2^ [2^. It seems that this state enjoys one more 



such distinction, this time in respect of our generalised 
EPR correlation. 

Conjecture 1 : Among all bipartite states of fixed entan- 
glement numerically equalling Er, and for every 9 in the 
range tanhr < tan 9 < 1, the two- mode-squeezed vac- 
uum I'i'r) yields the least value for the generalised EPR 
correlation Ag{-). In other words, no state {ip) with en- 
tanglement E{\ijj)) < Er can yield a generalised EPR 
correlation Ag{ip) < Ag{'i/r), for any 9 in the range 
tan 9 > tanh r 

The special case 9 = tt/A is the basis of the important 
work of Ref. [23|. Hence the present assertion can be 
viewed as a generalisation of their Proposition 1. 

The original EPR correlation A7r/4(-) continuously de- 
creases to zero with increasing entanglement. But this is 
not true of the generalised EPR correlation Ag{-). 

Let us denote by rg the value of r determined by a given 
value of 9 through the equation tan6' — tanhr, and let 
9r denote the value of 9 so determined by r. Then, for a 
given numerical Er, the relevant range for 9 in Conjec- 
ture 1 is 6*^ < 6* < 7r/4. 

Proposition 4: The generalised EPR correlation Ae(-) 
obeys the basic inequality Ae{-) > cos 26*. The two-mode- 
squeezed vacuum saturates this inequality if and only if 
the squeeze parameter r solves tanh r = tan 9. 
Proof: It is clear that the relations tan0 = tanhr, 
sin 29 = tanh 2r, and cos 29 = ( cosh 2r )^^ are equivalent 
to one another, and so also arc the inequalities tan^ > 
tanhr, sm29 > tanh2r, and cos2(? < (cosh2r)~^. Now 
consider the transformation (a, b) U{r){a, b)U{ry 
where U (r) — exp{ r (0^6^—06) } is the unitary two- mode- 
squeeze operation: 

a — > a cosh r — 6^ sinh r, b ^ b cosh r — sinh r. 

This implies the following transformation for the anti- 
commutator {b,b^ = bb^ -f 6^fe : 

_> {b'<b~ a^a) + ^{{a,a^ + {b,b^)cosh2r 

- {ab + a'' b'' ) sinh 2r 
= cosh2rAe^, 0^ = arctan(tanhr). 

Since {b,b^} > 1, so is also its unitary transform 
cosh2r Ae^. That is, Ag^ > (cosh2r)~^ = cos26'r. 

Thus, saturation of the inequality Ag^{ip') > cos29r 
is equivalent to the condition {4>\{b,b''}\tp) = 1, 
where = U{r)\tp). A pure state which satisfies 

{i>\{b,b'^}\^|:) = 1, is of the form IV') = \(f>) a (S) \0) b , where 
is any vector in Alice's Hilbert space Ha- It follows 
that states saturating the inequality Ag^{p) > cos29r 
constitute the set { p = U {r)pA ®\Q)b b (0|J7(r)^ }, where 
PA is any (pure or mixed) state of Alice's oscillator. Fi- 
nally, Conjecture 1 claims that among all these states sat- 
urating this inequality the two-mode-squeezed vacuum 
l^'rs), corresponding to the choice pA = |0)aa(0|, has 
the least entanglement. 
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Entanglement of Formation : With the canonical form 
and the generaUsed EPR correlations in hand, we are 
now fully equipped to compute the EOF of an arbitrary 
two-mode Gaussian state. 

Proposition 5 : Given an inseparable zero-mean two- 
mode Gaussian state pvb with covariance matrix Vq speci- 
fied in the canonical form by u, v, 9o and tq with u, v > 
and < tanhro < tan^o < 1, its EOF equals Er^, the 
entanglement of the squeezed vacuum \^ro)- 
Proof: The fact that M = Vq — V*^^ > guaranties 
that pvo can be realized as a convex sum of displaced 
versions D(^)|^ro) of the squeezed vacuum state \^ro): 
all of which have the same entanglement E^g as l^'ro): 



PVo 



Here D{^) is the unitary phase space displacement oper- 
ator. The rank of M equals 2, and both AI~^ and the 
two-dimensional integral refer to the restriction of the 
phase space variable ^ to the range of M. 

Since a specific ensemble realization with average en- 
tanglement Erg is exhibited, EOF(pvb) < Erg. On the 
other hand, evaluation of the generalised EPR correla- 
tion Ag{pvo) = tT(Aopvo), for the particular value of 
6 occurring in Vb shows that A^q (pvo ) = cosh 2ro — 
sin2^osinh2ro. And by Conjecture 1, this implies 
EOF(pvo) > Erg. We have thus proved EOF(pvi,) = Erg. 

An attractive feature of the canonical form of the co- 
variance matrix is that the two-mode-squeezing U (r) acts 
on it in a covariant or form-preserving manner. 
Proposition 6 : Under the two-mode-squeezing transfor- 
mation U (r) we have 



Vo{rQ,9o,u,v) 
r'o = rQ+ r, 
{u, v') 



Va{r'Q,e'^,u',v') ■ 

■ sinh2r 
sin 29'n = 



cosh 2r sin 2do 



cosh 2r + sin 29o sinh 2r ' 
{u, v) X (cosh 2r + sin 2^0 sinh 2r). 



This is easily verified by direct computation. While the 
canonical squeeze parameter rg simply gets translated 
by r, the parameters u and v get scaled by a common 
factor. If we define rgg , rg/^ through tan = tanh rgg and 
tan 9'q = tanh r^^ , the transformation law for takes the 
form of translation: r^^ = rgg + r. 

As a consequence of this covariance, the convex decom- 
position which minimizes the average entanglement goes 
covariantly to such a decomposition under two-mode- 
squeezing: the minimal decomposition commutes with 
squeezing. This implies, in particular, the following sim- 
ple behaviour of EOF under squeezing: Erg —^ Erg+r. 

Finally, the just separable Gaussian states on 
the separable-inseparable boundary, correspond to the 
canonical form with ro = As was to be expected, 
the condition (5) places no restriction on 9q in this case. 
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